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Traditional close-coupling methods suppose an expansion of the total wave func- 
tion in terms of inner stationary states of colliding subsystems. In the case of 
hadronic atoms, a similar expansion has to involve, inter alia, low angular momentum 
states (ns, np) with large annihilation or nuclear absorbtion widths. The life times 
T ni = h/^ni of these states are small as compared with the collision time and mean 
time between subsequent collisions, therefore the close-coupling approach has to be 
modified for the similar problems. 

In this paper we propose a generalization of the close-coupling method with anni- 
hilating states included in the basis. The correct asymptotic behaviour of the wave 
function in the annihilating channels suppose that the annihilating states can not be 
presented in the incoming channels whereas the corresponding components of the wave 
function of relative motion in outgoing channels have to damp out at large distances. 
The S-matrix of the transitions in the subspace a of all other (non-annihilating) states 
is not unitary, because the hamiltonian of the problem is non-hermitian. The unitary 
defect (1 — |<Sj-i| 2 ) gives the cross section of induced annihilation for the initial state 
i G a. The scheme of the numerical solving of the generalized close-coupling equations 
is proposed. It includes the calculations of two types of solution, from the origin to 
an intermediate sewing point R = a and from a large R to the point a. 



1 Introduction 

Quantum close-coupling method is used in the many theoretical problems of atomic and 
nuclear two-body collisions A + B — > A' + B' . It supposes an expansion of the total 
wave function ^(R, £) of the system in terms of inner wave functions of colliding 

subsystems, 



where the coefficients ^(R) have a sense of wave functions of the relative motion of two 
subsystems (A and B) in the corresponding channels. They satisfy to the system of cou- 
pled equation with the standard boundary conditions at R — > oo (incoming and outgoing 
waves in the open channels, and damping waves in the closed channels). The expansion (OQ) 
and the mentioned boundary conditions suppose that the basis states are stationary. 

However excited states of real quantum systems (atoms, nuclei, etc.) have finite life times. 
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Therefore the mentioned approach can be used if the life times of the states involved in the 
consideration are large as compared with the collision time, r k ^> r co n ~ Ro/v, where R is 
an interaction radius and v is the relative velocity of the colliding subsystems. This condi- 
tion, as a rule, is even not mentioned in the papers that used the close-coupling methods, 
because it is fulfilled in the most part of the problems of atomic and nuclear collisions. How- 
ever, there is, at least, one class of the problems where the mentioned condition is strongly 
violated, namely, when one of the colliding subsystems is an hadronic atom. In this case 
the expansion ([I]) has to involve, inter alia, low angular-momentum states (ns, np) with 
large annihilation (nuclear absorbtion) widths. The life times r n i = h/T n i of these states are 
small as compared with the collision time and mean time between subsequent collisions, and 
therefore these states can not be presented among possible incoming channels. They can 
be admixed to the total wave function in the area of interaction (R < R ), but then they 
disappear during collision. The corresponding components ^(R) of the wave function of 
relative motion in outgoing channels have to damp out at large distances. The parameter 

"nl = Tcoll/ T n/ ~ T n iR /tiV 

gives a criterion of the state stationarity during collisions. At v n \ <C 1 the non-stationarity 
during collision can be neglected, whereas for the states with v n \ > 1 the non-stationarity 
should be taken into account. We will refer similar short-lived states as annihilating states. 

The problems of very short-lived states in collisions are not specific for the close- coupling 
method. They would arise in any theoretical approach to the collisions of the similar 
systems. However they are considered up to now only in the framework of semiclassical 
approximations (see, e.g., [H [2]) that hold for "hot" hadronic atoms (ka ^> 1, or, for 
(pp)ni + H collisions, E > 1 eV). To our knowledge the single attempt to take into account 
nuclear absorbtion in rzs-states during collisions within the quantum close-coupling method 
was done in the paper [3] for the {^~p) n i + H scattering. However the authors use an 
artificial assumption that the widths T ns are turned off at some distance between two 
atoms (R> Ro = 5a ) that, evidently, contradicts to the physical reality. 

In this paper we propose a generalization of the close-coupling method with short-lived 
(annihilating) states included in the basis. The boundary conditions at R — > oo in the 
channels related to the annihilating states are formulated so that they have the correct 
quantum mechanical meaning in accordance with the physical properties of these states. 

2 Close-coupling equations with short-lived states in- 
cluded in the basis 

Let us consider multichannel two-body collisions of hadronic atom or ion (hZ) n i with a 
neutral atom A 

(hZ) nl + A^ {hZ) n , v + A, (2) 

where h is a negative hadron (ti~ , K~, p, etc.), Z is a bar nucleus (e.g., H + or He 2+ ). For 
the sake of simplicity we suppose that (i) the hadronic atom does not contain electrons, 
therefore the states of the (hZ) are classified by the principal quantum number n and the 
quantum numbers of orbital angular momentum I and m, (ii) the energy E of the collision is 
small (E <C E a , where E a = 27.2 eV is the atomic unit of energy) so that the channels with 
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the excited atom A* are closed and can be neglected, (iii) the principal quantum number n is 
large enough, therefore the orbital velocity of the hadron v n = Zv a /n <C v e , where v a = e 2 /h 
is the atomic unit of velocity and v e ~ Z e fjv a is a characteristic electron velocity in the 
atom A. The conditions on E and n mean that the heavy particles (nuclei and hadron) in 
the problem under consideration are slow as compared with electrons. Therefore electronic 
variables can be separated out within adiabatic approximation. Then the sum of Coulomb 
interactions between the hadronic atom (hZ) and the many-electronic atom A is replaced 
by the effective three-body interaction V(R, r) that depends on the relative coordinates 
R of two colliding subsystems and on the inner coordinates r of the hadronic atom. Total 
effective hamiltonian of the three heavy particles (h — Z — A) can be presented as 

H = T(R)+h(r) + V(R, r), (3) 

where T(R) = (— l/2m)VR is the kinetic energy operator, m is the reduced mass of the 
colliding subsystems, and h(r) is the inner hamiltonian of the hadronic atom. 

The hamiltonian of the hadronic atom could be, in general, very complicated due to 
possible transmutations of the (h + Z) system into other hadronic systems. However for 
our aims it can be presented as 

h = h + U apt (r), (4) 

where ho is the hamiltonian of hydrogen-like atom with the reduced mass \x and with the 
nuclear charge Z, and U opt (r) is a short-range complex optical potential with the absorptive 
imaginary part lmU opt (r) < 0. Let E ni and 0„; m (r) = R n i(r)Yi m (Q r ) be the eigenvalues 
and the eigenfunctions of the hamiltonian (111) , corresponding to discrete levels of the (hZ) 
system, 

h(p n i m (r) = E ni (f) nlrn (r). (5) 

The hamiltonian h is non-hermitian (h + ^ h) due to the imaginary part of the optical 
potential. Therefore the eigenvalue E n i in general case also can contain a negative imaginary 
part, 

E nl = E* - ilV2, (6) 

where T n i is the annihilation or nuclear absorption width of the level. The real part of the 
eigenvalue is usually expressed as E^ = e n — e n i, where e n = —fiZ 2 /2n 2 is the energy of the 
hydrogen-like system with the point nuclear charge, and (— e n {) is the energy shift of the 
level due to the nuclear interaction and the finite distribution of the nuclear charge. (We use 
the atomic system of units h — e — m e = 1, except as otherwise indicated explicitly.) Total 
difference (E n i — e n ) is usually referred as the complex energy shift AE n i = —e n i — iT n i/2. 

The conjugated hamiltonian h + has the complex-conjugated eigenvalues E n i = E* nl and 
the eigenfunctions n / m (r) = R n i(r)Yi m (fl r ) , where R n i(r) coincides with i?* ; (r). Two sets 
of the eigenfunctions, 0„z m (r) and n / m (r), form jointly a bi-orthogonal system, 

((pnlml^n'l'm') = {<t>nlm\4>n'Vm') = 8 nn >5u>5 mm r . (7) 

In order to introduce a close- coupling expansion of the total wave function of the system 
(hZ) + A, we define the basis states 

$/ M7r (r, = (0 ni (r) ® Y L (n R )) JM . (8) 
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Here Qr = (8, <fi) are polar and azimuthal angles of the vector R, L is the quantum number 
of the relative angular momentum of the colliding subsystems, and (<j> n i ® Yl)jm means 
the vector coupling of the states of two angular momenta (I and L) into the total angular 
momentum J. Index j = {n, I, L} enumerates the basis states with the definite quantum 
numbers of the total angular momentum J, M and of the parity 7r = (— l) l+L . If </> n zm(r) ^ 
4>nim{ r ), we have to introduce also the states $j M7T (r,Q R ) that are differ from (jSJ) only by 
substituting 4>ni{ r ) f° r 0n;( r )- I n this case we obtain again the bi-orthogonal system, 

^jm^j'm'^ = pswfijM^ = 5 l3 5 jr 5 MMI 5^, (9) 

Total wave function of the (hZ) + A system in the state with the total energy E, the total 
angular momentum J and the parity it in the framework of the close-coupling approximation 
can be written as 

*f JM7r (R, r) = J2 ^(r, n R )il£ J *{R)/R, (10) 

k 

where i is the index of an incoming channel. The conserved quantum numbers (E, J, M, it) 
will be omitted in the further equations the for the sake of simplicity 

Substituting the expansion fflQj) into the stationary Schroedinger equation with the 
effective hamiltonian ([3]) and taking into account the bi-orthogonality property (jSj), we get 
the system of coupled-channel equations 

^(R) + [k] - Lj(Lj + 1)/R 2 } 4j jt (R) = 2mJ2v jk (R)Ai(R), (H) 

k 

where fc| = 2m(E - E nl ) and V jk (R) = (<&j\V(~R, r)|$ fc ). The system (frTJ looks very 
similar to that one commonly used in the standard close-coupling methods. However in the 
problem under consideration the wave numbers kj in some channels can contain imaginary 
parts, therefore the conventional boundary conditions at R — ► oo (incoming + outgoing 
waves) in these channels have no sense and should be modified. For simplicity, we will not 
consider the closed channels with E^ > E. As for the behaviour of ipji(R) at R — > 0, the 
standard requirement of the regular solutions and the boundary conditions at the origin 
tpji(0) = remain valid. 

Let us divide the total space of the channels included in the basis into two subspaces, a 
and (3. The subspace a contains the channels that correspond to the stationary states with 
the negligible widths (T n i —>■ 0) and, respectively, to the very small values of the parameter 
u n i, whereas the subspace (3 contains annihilating channels {y n i > 1). As a rule, ns-states 
of hadronic atoms fall into the subspace (3. In addition, depending on the physical system 
and the kinetic energy under consideration, other states with small angular momenta can 
also fall into this subspace. The absorption widths T n i decrease drastically with angular 
momentum, therefore the states with large / always belong to the subspace a. The wave 
numbers kj for the channels j £ a are real, whereas for the channels j 6 (3 they are complex. 
In the latter case we define the wave number so that Imkj ^ 0. 

It is evident from the physical considerations that the channels from the subspace (3 can 
not be presented among the incoming channels, because the hadronic system (hZ) in the 
short-lived state will disappear before the arrival from a far distance into the interaction 
area. Therefore the corresponding components of the total wave function have to be zero 
everywhere in R G [0, oo): 

ij)ji(R) = at i G (3 at any j. (12) 
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(Recall that the second index in ipji(R) enumerates the incoming channels.) 

The remaining nonzero solutions ipji(R) correspond to the ordinary incoming channels 
(i £ a). The components with j G a can contain the both incoming (at j = i) and outgoing 
waves, whereas the components with j G ft have to damp out at the large distance. These 
latter can be presented in the form of 'outgoing' waves with the complex kj (Imkj ^ 0). 
Hence the asymptotic boundary conditions at R — ► oo for the functions ipji(R) at i G a 
and any j can be written as 

2 /<2 1/2 

ipji(R) — > h i exp[—i(kiR — Liir/2)]8ij — kj exp[i(k jR — Ljii/2)}Cji at i G a, (13) 

For the channels j G a this asymptotic expression has the traditional sense, and the co- 
efficients at the outgoing waves are the elements of S'-matrix, Cji = Sji, corresponding to 
the transitions between the channels in the subspace a. For the channels j G ft the first 
term in f|T3|) is absent due to j ^ i, whereas the second term ('outgoing' wave) is, in fact, 
damped out: 

ipji(R) ~ exp(— Imkj ■ R)Cji at % G a, j G ft. (14) 

The density current of the outgoing waves in these channels is also damped out exponen- 
tially, therefore the coefficients Cji at j G ft, i G a can not be treated as the S'-matrix 
elements of the transitions (a — > ft). 

It should be noted that the equation (1131) is valid in the far asymptotic region, where 
the centrifugal terms L(L + 1)/R 2 in the equations fill I) can be neglected. Instead, we can 
introduce the solutions in terms of the spherical Riccati-Hankel functions 

^ji(R) -> kT^h^ikiR^ij - kJ^h^ik^Cji (i G a, any j). (15) 

The spherical Riccati-Hankel functions are defined so that 

hf\x) -> exp[±i(x - Lrr/2)] 

at x — >• oo, therefore (JTBT) at very large R goes into ffl~3l) . But the solutions ffT5l) are valid 
also at the intermediate R outside the interaction well. 

Consider a possible scheme of the numerical solution of the system ffTTl) satisfying to the 
boundary conditions ipji(0) = at the origin as well as to the constraint f|T2|) at i G ft and 
to the boundary conditions (TLB"]) at R — > oo. Let iV be the total number of the channels 
included in the expansion ffTUl) . iV Q and Afg the numbers of the channels in the subspaces a 
and ft, respectively. The total number of the equations in the system ( TTT|) is also equal to N. 
Therefore the system of the second-order differential equations has iV linearly independent 
regular solutions f p (R) with the components f qp (R), where the first index q is the number 
of the component, and the second index p is the number of the regular solution. In the 
practical calculations the different regular solutions of the system are obtained, e.g., by 
using iV linearly independent initial values of the first derivatives fp(0). Total set of the 
regular solutions can be presented as the (N x N) matrix F(R) = {fi(R), f~2(-R), • • • , f/v(-R)}. 

Formally we could obtain the numerical solution for F(R) from the origin to the large 
R and then find the linear combinations of f p (-R) satisfying to the conditions f|T2l) and ([TBI . 
However this procedure leads to the essential difficulties. The imaginary part of energy in 
the annihilation channels produces the exponentially growing components of f p (R), which 
will be canceled in the correct linear combination of f p (R). But the remaining parts of the 
functions can not be calculated with a good accuracy at R — > oo. 
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In order to avoid a loss of accuracy, we use the solution F(R) only in the some restricted 
area R < a. For the area R > a we construct the solutions of the system fill I) beginning 
from a some external point R e > Rq, where the interaction Vji(R) is negligible, and going to 
the smaller R up to R = a. The choice of the sewing point a depends on the parameters of 
the physical problem under consideration. It has to be the subject of the numerical search 
in order to obtain a good accuracy of the results. 

Let yk{R) be the solution of the system ffTTl) that contains at a large distance the single 
outgoing or damped out wave in the channel k, so that 

vAfy^SjkkjWhfftoR). (16) 

Total set of the JV linearly independent solutions forms the (JV x JV) matrix Y(R) that 
tends to the diagonal matrix at a large R, but contains also non-diagonal elements in the 
interaction area (R < Rq). Similarly, we introduce the N a solutions Xj(i2) that contain at 
a large R the incoming waves in the channels % G a, 

Xji (R) ^Sjik-^h^ihR). (17) 

In line with the condition ffi~2l . this set of the solutions doesn't include the rising solutions 
at i £ /3. We introduce the square (JV x JV) matrix 



X(R) 



X aa 
Xp a 



that contains the N a columns of the Xj(i?) and the Np columns of zeros. 

Total (JV x JV) matrix of the linearly independent solutions of the system f ill I) satisfying 
to the condition ffT21) and to the asymptotic boundary conditions ([T2]) can be presented in 
the form 

V(R) = X(R)-Y(R)C, (19) 

The matrix ^(R) and the derivative ty'(R) have to be sewed at a point R = a with a linear 
combination F(R)A of the regular solutions, 

F(a)A = X(a) - Y{a)C, (20) 
F\a)A = X\a) - Y\a)C. (21) 

These equations allow to obtain the constant matrices A and C. Let P(R) and Q(R) be 
the (N x JV) matrix solutions of the system (TIT]) . Introduce the matrix Wronskian 

W(P T ,Q) = P T {R)Q'{R)-P IT {R)Q{R), 

where the superscript T means the transposition of the matrix. Due to the symmetry of 
the interaction in ffTT]) . Vjk(R) = Vkj(R), the Wronskian is a constant matrix, which is 
independent on R. For the solutions X(R) and Y(R) we obviously have 



W(X T , X) = W(Y T , Y) = 0, (22) 

la N 

0, 



W{Y 1 \X) = -2i( I « °Y (23) 
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where is the unit matrix in the subspace a. Using the equations (1201) and (I2T1) . we get 

A = W~ 1 (Y T , F)W(Y T ,X), (24) 
C = Y-\X-FA)= (^T [J), (25) 

where all solutions (X, Y, F) are taken in the sewing point R = a. 

The submatrix C aa , which gives the S-matrix of the transitions in the subspace a, can 
be expressed in terms of the Wronskian matrices: 

% = Cji = [W{X T ,F)W- 1 {Y T ,F)] j . at i,j e a. (26) 

The S-matrix is not unitary, because the hamiltonian of the problem is non-hermitian. The 
diagonal part of the unitary defect gives the probability of the induced annihilation for the 
incoming channel i6«, 

P(i -> annih) = 1 - {S^ 2 . (27) 

It allows to obtain the cross section of induced annihilation 

a? = (ir/k 2 )P(i — > annih). 

Total cross section of the induced annihilation is a sum of af with the relevant statistical 
weights of the channels. 



3 Conclusion 

We have generalized the close-coupling method with annihilating states included in the 
basis. The correct boundary conditions at infinity suppose that the annihilating states are 
absent among the incoming channels, but are presented among the outgoing channels as the 
damped out waves. The S-matrix of the transitions in the subspace a of all non-annihilating 
states is not unitary due to the non-hermitian hamiltonian of the problem. The unitary 
defect (1 — \Sji\ 2 ) gives the cross section of induced annihilation. 

The scheme of the numerical solution of the system of close-coupling equations is pro- 
posed. It contains the calculations of two types of solution, from the origin to a point R = a 
and from a large R to the point a, and the sewing of these solutions at R = a. This scheme 
gives the S-matrix of the transition between the ordinary channels as well as the probability 
of the induced annihilation. The choice of the sewing point a depends on the parameters 
of the physical problem under consideration. 

Application of the developed approach to the problem of collisions between antiprotonic 
helium ion (pHe 2+ )* and ordinary He atom at very low energy (~ 10 Kelvin) is given in 
the separate paper [3]. 
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